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We study the hopping conductivity in a composite made of Gaussian coils of flexible metallic 
wires randomly and isotropically suspended in an insulator at such concentrations that the spheres 
containing each wire overlap. Uncontrolled donors and acceptors in the insulator lead to random 
charging of wires and, hence, to a finite density of states (DOS) at the Fermi level. Then the Coulomb 
interactions between electrons of distant wires result in the soft Coulomb gap. At low temperatures 
the conductivity is due to variable range hopping of electrons between wires and obeys the Efros- 
Shklovskii (ES) law In cr cx —(Tes/TY^^ with Tes decreasing with concentration and length of the 
wires. Due to enhanced screening of Coulomb potentials, at large enough wire concentrations and 
sufficiently high temperatures, the ES law is replaced by the Mott law Intj oc —{Tm/T)^^'^. 



I. INTRODUCTION 

The electronic transport properties of composites made 
of metallic granules surrounded by some kind of insulator 
is widely studied. Spheroidal metallic particles or quan- 
tum dotsi^, relatively short carbon nanotubes^i^iSii, 
and long flexible nanowires such as bundles of carbon 
nanotubes or conducting polymersSiifliii are examples of 
the dispersed granules. If the concentration n of the 
metallic granules is large, they touch each other and the 
conductivity is metallic. When n is smaller than the 
percolation threshold ric, the composites are in the insu- 
lating side of the metal-insulator transition and hopping 
is the only mechanism of conduction. 

If each granule is neutral, adding or removing one elec- 
tron from it costs a considerable electrostatic energy. The 
charge transport is blocked by this Coulomb charging en- 
ergy and has activation dependence of the conductivity 
resembling the Arrhenius law. However it is generally 
observed that at low temperatures, the conductivity a 
exhibits a temperature dependent characteristic of vari- 
able range hopping (VRH) transport. 



<j = aocxp[-{To/T)% 



(1) 



where Tq is a characteristic temperature. For spheroidal 
particles or quantum dots^'^''^, a is observed to be 1/2. 
For bundles of single wall carbon nanotubes^, both a — 
1/2 (Ref. ii) and a = 1/4 (Ref. ii^) are reported. 
In composites made of conducting polymers dissolved in 
some kind of insulating polymer matrijiSiiSiii, the index 
a was found to be 1/2 below the percolation threshold 
and at larger concentrations, a evolves from 1/4 to 1 
(Ref.^). The finite DOS at the Fermi level necessary for 
the VRH conductivity is attributed to the charged im- 
purities in the insulator playing the role of randomly bi- 
ased gates. Then the interaction of electrons residing on 
different wires creates the Coulomb gap, which leads to 
ES VRH between granulesii^. In our recent paper—, we 
studied the puzzling crossover from ES law with a = 1/2 
to Mott law with a = 1/4 for a suspension of straight 
metallic wires. We showed that the characteristic tem- 
perature Tes oc l/{nL^)'^, where L is the wire length. 



In this paper we would like to calculate the hopping 
conductivity for the more complicated case of flexible 
nanowires suspended in the insulating medium. They 
can be just very thin metallic wires suspended in an in- 
sulating liquid which was abruptly frozen, or very long 
single wall nanotubes or their bundles, or heavily doped 
conducting polymers. Such objects can have one or many 
parallel metallic channels. In the latter case, particularly 
for doped conducting polymers, disorder should be im- 
portant and leads to localization of states at some dis- 
tance along the nanowire^^-. Nevertheless in this paper 
we concentrate on a fundamental problem of clean or 
weakly disordered wires. Generalization of our theory 
to the case of strongly disordered wires, where electrons 
are localized within distances much smaller than the wire 
length, will be discussed in the conclusion. 




FIG. 1: Local view of the semidilute solutions of flexible Gaus- 
sian wires. We show the case where the mesh size is not longer 
than persistence length, the wire within each mesh is essen- 
tially straight. At lesser density, the mesh size is longer than 
the persistence length, and then the wire in the mesh is wiggly 
(not shown). 

We assume that each wire has the contour length L, 
relatively large persistence length p and the diameter 
(thickness) d. We assume L ^ p ^ d. Throughout 
this work we require L < p^/d^ and this lets us disregard 
the effect of excluded volume, considering each wire as 
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a Gaussian coil with the size R = y/Lp. At small con- 
centrations, where the inter-coil distance is much larger 
than the coil size, we can simply utilize the results for 
arrays of quantum dots - . The overlap of the coils starts 
when the concentration of wires exceeds the concentra- 
tion of the wire inside one coil: n — l/{Lp)^^^ (see Fig. 

Therefore our discussion is confined to the range of 
concentrations l/{Lp)^^'^ ^ n ^ ric, where Uc = l/L^d 
is the percolation threshold at which the hoppping con- 
ductivity is replaced by the metallic one (we derive the 
expression of below). 

In this paper we restrict ourselves to the scaling ap- 
proximation for the conductivity and to delineating the 
corresponding scaling regimes. In our scaling theory, we 
drop away both all numerical factors and, moreover, also 
all logarithmic factors, which do exist in the problem, be- 
cause it deals with strongly elongated cylinders. In this 
context, we will use symbol "=" to mean "equal up to a 
numerical coefficient or a logarithmic factor" . 

Our main results for the hopping conduction are pre- 
sented by the phase diagrams in the plane of parameters 
n(Lp)3/2 g^j^j y shown in Figs. EElandEl Fig. repre- 
sents the simplest case of narrow wires with only one 
conducting channel and relatively long length p^/d < 
L < p^/d^. Remarkably the ES hopping regime which 
vanishes for single channel straight wiresi^ does exist for 
Gaussian coiled wires. The reason is that each straight 
wire is a one dimensional object where Coulomb inter- 
action plays only marginal role. On the other hand, a 
coiled metallic wire from the point of view of the Coulomb 
interaction is a two-dimensional object, where charging 
energy Ec can be substantially larger than the quantum 
level spacing S. At large enough n{Lp)^^'^ and sufficiently 
high temperature, ES hopping conductivity is replaced 
by the Mott law due to increasing screening of coulomb 
potentials (regime M). At even higher temperatures both 
VRH regimes crossover to the nearest neighbor hopping 
(NNH) regimes. 




FIG. 2: (Color online) Summary of scaling regimes for the 
case of flexible Gaussian coiled wires each with single con- 
ducting channel and length p^/d < L < jS . Dark line sep- 
arates regimes where Coulomb interaction is important from 
regimes where it plays minor role. The grey (blue) line sep- 
arates activated NNH regimes Ni and N2 from VRH regimes 
ES and M. 



The plan of this article is as follows. In Sec. ^ we 
first consider the macroscopic dielectric constant Koff of 
the suspension. We then calculate the finite DOS at the 
Fermi level in Sec. IIIII We continue in Sec. IIVI bv cal- 
culating the effective localization length ^ for hop over 
distances larger than the typical coil size. Using these re- 
sults, we summarize all the scaling regimes for the hoping 
conductivities in Sec. Finally in Sec. IVII we conclude 
with the discussion of generalization of our theory to the 
case where wires are strongly disordered. 



II. DIELECTRIC CONSTANT 

We start with determining the static macroscopic di- 
electric constant Keff of the suspension with wire concen- 
tration Xjl?' <^ n ^ ric- The percolation threshold ric 
can be estimated from the following argument. Perco- 
lation happens when each wire has roughly two direct 
contacts with other wires. Following Ref.^^, the number 
of contacts per wire is of the order of nL'^d (each piece 
of the wire of the length d has probability nLd^ to be 
touched). Requiring that nL^d is of the order of 1, we 
obtain the percolation threshold ric — 1/L^dJ^ 

In Refii^, we have calculated the macroscopic conduc- 
tivity cr of a suspension of Gaussian coiled wires with con- 
ductivity ai in a weakly conducting medium with con- 
ductivity (72 <C cTi. When dealing with metal wires in 
an insulator, one can consider small final frequency re- 
sponse and use — icon/ At: for (T2- This lets us find the 
macroscopic effective conductivity a{uj) and Kes using 

KcS — 4:TTia{LL!)/uj, 

, (2) 

nnL^p it ^< n{Lpfl^ < 

where k is the dielectric constant of the insulator. We 
notice that the polarization of these metallic wires leads 
to great enhancement of the macroscopic dielectric con- 
stant Koff much larger than n of the insulating host. The 
derivation of Eq. ||2Jl is based upon the geometrical prop- 
erties of the system, therefore let us dwell on them first. 

In the parlance of polymer science, when wires strongly 
overlap at n > l/lLp)^^"^, we deal with a semidilute solu- 
tion - the system locally looks like a network with certain 
mesh size r (see Fig. In the scaling sense, r is the 
same as the characteristic radius of density-density cor- 
relation function, and can be estimated following Ref^. 
Suppose that the wire within each mesh has a contour 
length g. It makes a density about g/ Lr^ which must be 
about the overall average density n. Thus, g/Lr^ ~ n. 
Second relation between g and r depends on whether 
mesh size is bigger or smaller than persistence length p: 

i y/gp 9> P 
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Accordingly, one obtains 

9 = n^L^pa , r — ^ji^ if jj^^ < n < j-^ 
9 = \^ , r = if 1 < n < 4;. ' 

The upper line corresponds to such a dense network that 
every mesh is shorter than persistence length and wire is 
essentially straight within each mesh (see Fig. The 
lower line describes much less concentrated network, in 
which every mesh is represented by a little Gaussian coil. 

Now we present a simple derivation of Eqs. (jSJ equiva- 
lent to what we have done for the macroscopic conductiv- 
ity in Ref.^^. Consider a cube of the size about (Lp)^/^ 
inside our macroscopic sample. On the one hand, the 
overall dielectric constant on the scale of this cube is al- 
ready about the same as that of a macroscopic body, we 
denote it Kcff- Therefore the capacitance of such a cube 
is about (Lp)^/^Kcff . On the other hand, we can estimate 
this capacitance considering wires inside the cube. Be- 
cause the distance between wires is about r, each wire 
bridges another wire through a blob of insulator with 
size r and capacitance tk. Since the contour length of 
the wire inside each blob is about g, there are about L/ g 
such capacitors connected in parallel, yielding overall ca- 
pacitance connecting the given wire as {L/g){rK). Since 
all (or sizeable fraction of all) n{Lp)^^'^ wires in the cube 
are in parallel, we finally arrive at the cube capacitance 
as n{Lp)^/^{L/g){rK). Equating this to (Lp)^/^Koff and 
accounting for the different geometrical properties repre- 
sented by Eqs. 10}, we arrive at Eq. ||2J). 

Eqs. can also be understood along the following 
way. If the wave vector q is larger than l/{Lpy^^, the 
static dielectric function for the wire suspension has the 
metal-like form 

e(g) = K (^1 + , (5) 

where r represents the typical separation of the wire 
from other wires and plays the role of screening radius 
for the wire charge. The function e{q) grows with de- 
creasing q until q = l/{Lpy^^ where the composite loses 
its metallic response and e{q) saturates. As a result, 
the macroscopic effective dielectric constant is given by 
Koff = ({q = i/iLpY^^). Plugging proper r at different 
concentrations (Eq. Q), we obtain Kofi represented by 
Eqs. (121). They show how Kcs grows with L at L ^ p for 
flexible nanowires. In the following sections, we will use 
Eqs. Q for the calculation of the hopping conductivity. 
Note, however, that the dielectric constant of our system 
may be itself of interest in applications where very large 
dielectric constants are welcomed. 



III. DENSITY OF STATES 

As we mentioned in the introduction, according to 
ReffS, the finite density of states (DOS) near the Fermi 



level originates from uncontrollable donors (or acceptors) 
in the insulating host. Donors have the electron energy 
above the Fermi energy of wires. Therefore, they donate 
electrons to wires. A positively charged ionized donor can 
attract and effectively bind fractional negative charges 
on all neighboring wires, leaving the rest of each wire 
fractionally charged. At a large enough average num- 
ber (3> 1) of donors per wire, effective fractional charges 
on different wires are uniformly distributed from — e/2 to 
e/2. In such a way the Coulomb blockade in a single wire 
is lifted and the discrete density of states get smeared. 
If the charging energy Ec for a coil is larger than the 
mean quantum level spacing S inside the wire, the DOS 
go is given by 50 = ^/EdLp)^/"^ . Substituting the effec- 
tive dielectric constants given by Eqs. Q, we obtain the 
charging energy 




eV(Kn'i'/V/') if 1 < n(Lp)3/2 < ^ 
ey{nnL^/'py') if ^ < n{Lpf/^ < 



In the very vicinity of the Fermi energy, the long range 
Coulomb interaction creates the parabolic Coulomb gap 
A which leads to the ES VRH. 

If Ec is smaller than the mean quantum level spacing 
S, the Coulomb interaction plays minor role. The density 
of states go becomes l/S{LpY^'^. Level spacing S strongly 
depends on the cross section of the wire. If the wire is 
so narrow that it has only one conducting channel, 6 is 
the order of /mXpL, where Xp is the Fermi wavelength 
and m is the effective electron mass. In order to compare 
with the charging energy E^ we introduce constant P = 
Xp/a and rewrite 6 as c'^/PkL. Most likely, (3 is the 
order of 1 . If the wire is thick enough, the electrons have 
three dimensional energy dispersion and S is the order 
of Xph^ /rnL(P or Pe^a?' / kcPL. For both cases, the ratio 
Ec/5 decreases with growing n or i because of enhanced 
screening of Coulomb interaction in a denser system of 
longer wires. To calculate the VRH conductivity, we need 
now another key quantity-the localization length ^. 

IV. LOCALIZATION LENGTH 

It was suggestedii that long range hopping process ex- 
ceeding the size of a single wire can be realized by tun- 
neling through a sequence of wires. The states of the 
intermediate wires participate in the tunneling process 
as virtual states. The virtual electron tunneling through 
a single granule is called co-tunnelingy^iMMi and re- 
garded as a key mechanism of low temperature charge 
transport via quantum dots. One should distinguish the 
two co-tunneling mechanisms, elastic and inelastic. Dur- 
ing the process of elastic co-tunneling, the electron tun- 
neling through an intermediate virtual state in the gran- 
ule leaves the granule with the same energy as its initial 
state. On the contrary, the tunneling electron in the in- 
elastic co-tunneling mechanism leaves the granule with 
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an excited electron-hole pair behind it. Which mecha- 
nism dominates the transport depends on the temper- 
ature. Inelastic co-tunneling dominates at T > ^/Ec5, 
while below this temperature, elastic co-tunneling wins. 
In discussions below, we consider sufficiently low temper- 
atures such that only elastic co-tunneling is involved. 

Let us remind the idea of the calculation of the local- 
ization length ^ used in Ref&iLSS. When an electron 
tunnels through the insulator between wires at the near- 
est neighboring hopping (NNH) distance r„, it accumu- 
lates dimensionless action r„/a, where a is the tunnel- 
ing decay length in the insulator. We want to empha- 
size that the NNH distance r„ is realized only in one 
point of the wire and, therefore, is much shorter than 
the typical separation along the wire from other wires 
r. NNH distance r„ can be calculated using the percola- 
tion method^'^. If one consider each wire with thickness 
r„, using the same idea we showed above, percolation 
through these cylinders appears, when nL^r„ = 1. Thus 
r„ = 1/nL^. Since we assume the metallic wire is only 
weakly disordered, we neglect the decay of electron wave 
function in the wire. Over distance x, electron tunnels 
x/^LpY^"^ times accumulating a dimensionless action of 
the order of r„/a each time. Thus, its wave function 
decays as exp[— a;r„/(Lp)^/^a] — exp(— x/^), where the 
localization length 

^ = a(Lp)i/Vr„ = a(nL5/2pi/2). 

Now we have all the elements needed to construct our 
phase diagrams of the conductivity at various tempera- 
tures and concentrations. 



V. SUMMARY OF PHASE DIAGRAMS 




FIG. 3: (Color online) Summary of scaling regimes for the 
case of flexible Gaussian coiled wires each with single con- 
ducting channel and length p < L < / d. The dashed line 
marks the border where the wire within each mesh changes 
from Gaussian coiled to be straight. The efi'ective dielectric 
constant also changes at this line. 
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FIG. 4: (Color online) Summary of scaling regimes for the 
case of thick wires with length p < L < /d. 



We first consider narrow wires each with only one con- 
ducting channel. Depending on the relation between L 
and p, percolation can start before or after n{Lp)^^^ — 
(i/p)^/^, where the Gaussian coiled wire within each 
mesh becomes straight at higher concentrations. Our 
results for these two cases are summarized in Figs. |2 
and 121 in the plane of parameters n{Lp)^^'^ and T. Each 
line represents a smooth crossover except the percolation 
threshold n{LpY^^ = p^/^/L-^/^d, where the insulator- 
metal transition happens. In the vicinity of this transi- 
tion, the conductivity has a critical behavior, which goes 
beyond the scope of our scaling approach. 

To be systematic, let us start the review of scaling 
regimes from Fig. |21 which is valid when p^ /d < L < 
p'^ /d^ and therefore percolation happens at relatively di- 
lute solution with Gaussian coiled wire inside each mesh. 
We begin with activated NNH hopping regimes Ni and 
N2 at sufficiently high temperatures. The temperature 
dependence of the conductivity is given by 



where the first term is due to electron wave function 
overlap and in the second term Ea ~ max{Ec, 6} is the 
activation energy at very small tunneling conductances. 
Here the charging energy Ec should be taken from the 
upper line of Eq. JBJ. Due to enhanced screening of 
Coulomb potentials with growing n(Lp)'^/^, E^ decreases. 
As a result, the activation energy Ea evolves from Ec in 
regime Ni to the quantum level spacing 5 in regime N2 
at n{Lp)^/^ = ^i/2Li/4/pi/4^ j,^^ q jg ^^jj^ 

far as the contribution of the activation is smaller than 
that of the overlap. When they are comparable at suf- 
ficiently low temperatures, the NNH regimes smoothly 
crossover to VRH regimes at T = e^a/ nriL^^'^p'' ^'^ and 
T — e^anL/ (3k for the regimes Ni and N2 respectively. 

In the regime ES, Coulomb interaction is important, 
which leads to the Efros-Shklovskii (ES) VRH conduc- 
tivity 

a = aacM-[TEslTY'\ (9) 
where Tes = e^/KoffC- Applying the above line of Eqs. 
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(0) and Eq. lO, we obtain 



Tes 



l3Lll/2p7/2a 



(10) 



At large n{LpY/'^ and relatively high temperatures, the 
ES law is replaced by the conventional Mott law, 



(11) 



where Tm — 1 /goC^ ■ Plugging in go for the narrow wires 
and Eq. (O, we arrive at 



T, 



M 



(12) 



The regime ES borders the regime M along the line 
T — (ie^a/ Kn^L'^p'^ , which can be obtained by equating 
the Coulomb gap width A = \/goe^/f^ to the Mott's 
optimal band width T^^^/igoS.^y^'^- 

When the wire is so short that p < L < p^/d, perco- 
lation happens at relatively dense system where the wire 
within each mesh is straight. As a result, to the right of 
the transition from the coiled mesh to the straight mesh 
at n(Lp)'^/^ = {L/pY/"^, the effective dielectric constant 
should be taken from the bottom line of Eq. ((SJ. Thus 
for the regime ESi, the characteristic temperature Tes 
in Eq. should be replaced by 



TeSi = 



KarP'L'^l^p^l'^ 



(13) 



The crossover from this new regime ESi to the regime 
M is along the line T — fie^aj Knl?p^ . It is worthwhile 
to emphasize that Coulomb interaction loses its impor- 
tance for the NNH regimes at n{Lp)^^'^ = l3'^/'^{L/p)'^/^ 
which is smaller than (L/p)^/^, where the geometry of 
the mesh changes and thus so does the effective dielec- 
tric constant kcs- As a result, compared to Fig. [21 the 
regime ESi is the only new regime added in Fig. O At 
p = L, Fig. 13 resembles the case of straight wires, the 
ES regimes almost vanish. Until now, we complete the 
review of phase diagrams for the narrow wires each with 
single conducting channel. 

For thick wires, simply using the corresponding den- 
sity of states discussed in the Sec. IIIII we obtain phase 
diagrams similar to those for the case of narrow wires. 



but with new borders between scaling regimes. Unlike 
narrow wires, there is a possibility that the geometry of 
the mesh changes when the Coulomb interaction is still 
important. Scaling regimes of this case is shown in Fig. 
2] There is a new NNH regime N3 where the the conduc- 
tivity is represented by Eq. (jS)) with activation energy 
Ea given by the bottom line of Eq. ©. Replacing p 
by L in Fig. ^ we again recover the result we got for 
straight wires in Fig. 2 of RefiiS. 



VI. CONCLUSION 

In this paper we calculated the macroscopic dielectric 
constant and the hopping conductivity of a suspension 
of Gaussian coiled metallic wires in an insulator. We ob- 
tained a plethora of different scaling regimes. Our results 
are applicable to wires made of well conducting metals 
with relatively small disorder. They can also be general- 
ized to the case of strongly disordered wires where elec- 
trons are localized at a distance A much smaller than the 
length of the wire L. To this end, one should proceed sim- 
ilar to what we have done for straight wires in the Ref. «. 
To calculate K^ff, one can imagine that the disorder ef- 
fectively cuts the wire into shorter metallic pieces each 
with length A. Therefore, we can calculate KcS thinking 
about a suspension of metal wires with shorter lengthes 
A but larger concentration nL/X. Resulting KeS depends 
on A and is much smaller than that of clean wires. Calcu- 
lating the effective localization length ^ similarly to Sec. 
IIVI we should add the decay of electron wave functions 
along the wire. Taking into account the fractal proper- 
ties of the wire, we arrive at a much smaller localization 
length than that of clean metal. Reduction of Kcs and ^ 
leads to much larger Tes, Tm and much smaller hopping 
conductivity. We can construct phase diagrams which 
topologically look similar to diagrams of Figs. |2|21and 
^ but they are more complicated and we are not showing 
them here. 
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